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We study spatial correlations and structure factors in a three-state stochastic lattice gas, con- 
sisting of holes and two oppositely "charged" species of particles, subject to an "electric" field at 
zero total charge. The dynamics consists of two nearest-neighbor exchange processes, occuring on 
different times scales, namely, particle-hole and particle-particle exchanges. Using both, Langevin 
, equations and Monte Carlo simulations, we study the steady-state structure factors and correlation 

■ functions in the disordered phase, where density profiles are homogeneous. In contrast to equihb- 
' rium systems, the average structure factors here show a discontinuity singularity at the origin. The 

associated spatial correlation functions exhibit intricate crossovers between exponential decays and 
_i power laws of different kinds. The full probability distributions of the structure factors are universal 

I asymmetric exponential distributions. 

D 

B 

' I. INTRODUCTION 

The study of correlations and structure factors provides a sensitive probe into the characteristics of collective 
behavior in many-particle systems. For example, in a system with short range microscopic interactions, maintained 
' in thermal equilibrium, spatial correlations will in general decay exponentially. Long-range spatial correlations, 
characterized by power laws, are observed only if the system is at a critical point. In contrast, when such systems 
O ' are driven into non-equilibrium steady states, long-range correlations are often present in large regions of the phase 
^ O ^ , diagram ^ . A model system in which such anomalous correlations are easily observed and studied is the driven Ising 
lattice gas (the "standard model") [^,1). The drive ("external field") biases particle jumps along a specific lattice axis, 
fS| , thus forcing the system into a non-equilibrium steady state. One of the most intriguing and unexpected features of 
^ ■ this system is the presence of long-range spatial correlations at all temperatures above criticality, due to the breakdown 
[ of the traditional fluctuation-dissipation relations Q in conjunction with a conservation law for the particle density. 
] . In momentum space, the structure factor develops a discontinuity singularity at the origin 
_q ' In experiments, correlations are typically studied by photon, electron or neutron beam scattering techniques. The 
I , scattering intensity is closely related to the structure factor. Depending on the actual physical system, this quantity is 
— . ■ the Fourier transform of the "density-density" correlations, where, e.g., in a ferromagnetic system "density" stands for 
^\ , the local magnetization. Even in the stationary case, i.e., when the averages are not expected to be time-dependent, 

■ the densities themselves are fluctuating quantities in both space and time. Thus, when data are taken, it is crucial 
C3 ' to compare the time scale of these fluctuations to the duration of a typical "snap-shot" . If the former is much 
^ . shorter than the sampling time interval, then even one measurement practically results in a temporal average. In 

^ ] this case, the scattering intensity is a direct measure of the average structure factor. In the opposite scenario each 
individual snap shot appears as a random pattern of speckles. The collection of many snap shots, however, represents 
the full distribution of the fluctuating density-density products. This phenomenon has long been known in laser 
Q ] scattering experiments and the statistical properties of the random speckles are well established |^ . Using Monte Carlo 
simulations, it is particularly easy to probe fluctuating quantities in terms of the their distributions: each measurement 
corresponds to one configuration at a certain instant of time so that there are no "experimental" difficulties in 
achieving fine sampling. Since driven lattice gases exhibit generically singular density-density correlations, a study of 
5— i ' the associated distributions is expected to be particularly interesting. 

In this paper, we will focus on a generalization of the standard model, similar in spirit to the one leading from the 
Ising model to spin-1 |Q or Potts |^ models. Instead of just a single species of particles, we now consider two (labelled 
as -f 's and — 's) which are driven in opposite directions, subject to periodic boundary conditions. Empty lattice sites 
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are referred to as holes. This generalization is motivated by a variety of physical systems, ranging from fast ionic 
conductors with several mobile ion species |^ and water droplets in microemulsions with distinct charges ||ic|] to gel 
electrophoresis and traffic flow [Q. For simplicity, we neglect the usual Ising nearest neighbor interaction and 
retain "only" the excluded volume constraint. The model thus corresponds to the high-tempcrature, large-drive limit 
of a more complicated interacting system. 

This multi-species model, in both one and two dimensions, has been studied in detail In its simplest version, 

particles are allowed to exchange with holes only. Monte Carlo simulations [ p^ in two dimensions and mean-field 
studies |p^ , p^ show that there is a transition, controlled by particle density and drive, from a spatially homogeneous 
(disordered) phase to a charge segregated one, where the excluded volume constraint leads to the mutual blocking of 
particles. In this paper we extend our previous studies on structure factors and their distributions |Q to the 
more general case where the particles are also allowed to exchange amongst themselves: we "soften" the excluded 
volume constraint by allowing exchanges of nearest neighbor, oppositely charged particles on a time scale 7 which 
is distinct from that of the particle-hole exchange Here, the blocking transition still occurs, as part of a more 
complex phase diagram . A particularly interesting case emerges when the two time scales are chosen to be equal, 
i.e., 7=1: here, egzta/ charges are completely uncorrelated (up to trivial finite-size effects) while "hidden", non-trivial 
correlations survive between opposite charges. We add that this model in one dimension, at infinite drive and for 
arbitrary 7, has been solved exactly by Godreche and Sandow |p^ . 

We will focus on the disordered phase of the system, where we have a sound analytic understanding of the dynamics 
in terms of Langevin equations. We will study not only the averages but also the full distributions of the steady-state 
structure factors, using Monte Carlo simulations and a continuum field theory. Finding excellent agreement between 
our simulations and analytic results, we can trace the characteristics of the distributions back to the structure of the 
underlying Langevin equations. Given these relations, measurements of structure factor distributions in real systems 
should provide considerable information about the associated dynamics. 

The paper is organized as follows. In the next section, we define the microscopic model and give some details of 
the simulations. In Section III, we introduce the underlying Langevin equations and calculate the average structure 
factors, the corresponding spatial correlations, and the theoretical distributions of the structure factors. In the final 
section, we discuss our results and conclude with a brief summary. 



II. THE MICROSCOPIC MODEL 



We consider a two dimensional fully periodic lattice with L x L sites, each of which can be empty or occupied 
by a single particle. To account for the presence of two species, we introduce two occupation numbers and n~, 
with n being or 1, depending on whether a positive or negative particle is present at site x. The excluded volume 
constraint imphes n^n~ = , for any x. To model the system at zero total charge, we choose X^xt'^x ~ '^xl = 0' i-^-' 
the average densities of positive and negative particles are the same: 

P= ;^5I^x = ;^5I"x , (1) 

X X 

Since the dynamics conserves both densities separately, p is a constant. In the absence of the drive, the two species 
of particles are distinguished only by their label: both types hop randomly to nearest-neighbor empty sites, with the 
same rate F. In addition, nearest neighbor pairs of opposite charges exchange with a rate 7F. The external drive 
is directed along a specific lattice axis, labelled as the -fX||-direction. Reminiscent of a uniform "electric" field E, 
it exponentially suppresses jumps against the force. Specifically, during one Monte Carlo step 2L^ nearest-neighbor 
bonds are selected at random. If a particle-hole pair is encountered, an exchange takes place with probability 

Wph ^ T mm{l, exp{qE Sxi\)} , (2) 

where q = ±1 is the charge of the particle and (5x|| = ±1, is the change of the 2;|| coordinate of the particle due to 
the jump. Similarly, if the neighboring sites are occupied by opposite charges, a particle-particle exchange (or charge 
transfer) is attempted with probability 

Wpp = 7Fmin{l,exp(i;(5x||)} , (3) 

where now Sx\\ is the change in the coordinate of the positive particle due to the jump. Note that we do not 
introduce a factor of 2 in the exponential here, as one might have expected for a real electric field. This choice leads 
to a simpler Langevin equation without significantly affecting the phase diagram. Needless to say, it is irrelevant 
whether exchange takes place or not, if both sites carry identical content. 
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For our simulations, we set F = 1, so the control parameters are p, E and 7. On lattices with L ranging from 30 
to 100, the system is initialized with random configurations of various particle densities. Runs last from 2.5 x 10^ to 
5 X 10^ MCS. The first 62500 MCS are discarded to allow the system to settle into steady state. Then, we measure 
the Fourier transforms of every 125 MCS, defining them in the usual way: 

^E^^^'^'^x • (4) 

X 

In the following, we will investigate equal-time density-density operators in momentum space, considering both their 
full distributions as well as their (ensemble or time) averages. In the literature, the term "structure factor" typically 
refers to the averages, i.e., 

5"^(k)^iKn!j, (5) 

where a, /? = +,— ; k = ^(m^, TO||) 7^ and V = is the volume. Occasionally, especially when discussing the full 
distributions, we will use the word "structure factor" for the fluctuating two-point operator itself. In the disordered 
phase, S"^ is the Fourier transform of the usual equal-time correlation function 

G"''(x)^(n>f)-«)(nf). (6) 

Thus, if G is even in x, S will be real, so that an imaginary part of S signals a part of G which is odd in x. By charge 
symmetry, we expect G^^ = G . Clearly, both must be even in x, so that the associated S"s are real. On other 
hand, we have 

G+-{^,E)=G+-{x^,-xp-E) (7) 

in the presence of the drive so that may have an imaginary part (which must be odd in E) . Finally, G ^ (x) = 
G^ (^x) follows from (^) by translation invariance. Turning to the full distributions, these can be constructed from 

the time series of y , ^ — — and ^ — — m the steady state, ijxpfoitmg symmetries again, we note that 
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y and '^y ^ are distributed identically, so that we need to consider only the former. Further, only ''^ ^ and 

"''y"'' are necessarily real, while ^y''^ will generically be complex. 

We simulate systems with 7 ranging from to 1. For small 7's we choose E and the density p in such a way that the 
system is in the homogeneous phase. For larger 7's (7 > 7c ~ 0.62) the charge exchange mechanism suppresses the 
ordered phase entirely [p^ so we can pick arbitrarily large fields at any density. A particularly interesting case occurs 
for 7=1. Here, the rates for particle-hole and particle-particle exchanges become equal, i.e., Wpp = Wph, so that a 
positive (negative) particle can no longer distinguish a negative (positive) one from a hole. Thus, a positive (negative) 
particle experiences biased diffusion, slowed only by encounters with other positive (negative) particles, just as in the 
case of a single, non-interacting species, whose steady state probability distribution of configurations (i.e., the steady 
state solution of the associated master equation) is exactly known to be uniform (21] . For our case, this implies that the 
marginal distribution of the occupation numbers of one species is uniform, i.e., P [{n^}] = X]{nT} P [{n'^,n~}] oc 1. 

Thus, we expect questions regarding only one species of particles to have trivial answers, e.g. G^''~(x) must vanish 
for X 7^ in an infinite system or yield the finite-size fluctuations in a finite one. On the other hand, the two-point 
function between opposite charges can display interesting structures, e.g. long-range correlations, as a result of the full 
distribution P n^}] not being uniform. We note briefly that a completely "flat" steady state, P n~}] cx 1, 
is obtained for 7 = 2, as in the one-dimensional version of our model [ p8| . 

In Fig. |l|, we present the results for the three independent S"s found in the 100 x 100 system at a small value of 7 
and note the discontinuity singularity of these objects at the origin. In Fig. |[ we show the same quantities for 7 = 1 
and draw special attention to the fact that, while S'^'^ does not depend on k at all, exhibits a highly nontrivial 
k-dependence. Fig. |^ and ^ present the structure factor distributions for the smallest longitudinal and transverse 
wave vectors, respectively. Before discussing the data in detail, we will first present the theoretical framework within 
which they can be understood. In particular, we will focus on two points, namely first, the emergence of discontinuity 
singularities in the structure factors at k = 0, and their consequences for long-range correlations in real space, and 
second, the origin of the asymmetric exponential form of the distributions. This will then be followed by a comparison 
between our theoretical predictions and the simulations. 
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III. COARSE-GRAINED DESCRIPTION 



To extract the behavior at large distances (or small k in momentum space), a continuum field theory for the slow 
variables of the model is most appropriate. To find such a description, we must (i) identify the slow variables of the 
theory, and (ii) obtain a set of equations of motion for these quantities, corresponding to a coarse-grained version of 
the microscopic dynamics. For generality, we consider the d-dimensional case when X|| is directed along the electric 
field and xj^ is in the (d — l)-dimensional subspace, perpendicular to the field. Time is denoted by t. Then, the 
slow variables are easily identified as the conserved densities, p^(x, t), of the two species. The most systematic way 
to arrive at their equations of motion is to perform an Jl-expansion p2| , ^ : After partitioning the whole system 
into sufficiently large blocks of size il, one splits the particle densities associated with the block centered at x into a 
macroscopic part (p^) and a fluctuating one (x^) : 

^ E "± =p±(x,t) + f]-VV(x,t). (8) 

x' eb(x) 

This decomposition is inserted into the microscopic master equation, followed by a systematic expansion in fi. At 
leading order, we find a set of mean-field equations of motion for the local densities which reads, after taking a naive 
continuum limit: 

dtp^ = -Vr I [p± V (1 - p+ - p-) ± £X||p±(l -p+- p-)] + 7[p± V ± £X||p±pT]| , (9) 

where 



r 

r,. 



(10) 

II / 

is the diffusion-matrix. Tj_ is diagonal and isotropic in the {d — l)-dimensional subspace, thus characterized 

by a number Fj^. V is the asymmetric gradient operator, acting on any two functions / and g according to 

f 'V g ~ f^g — g^f- £ is the coarse-grained bias and X|| is the unit vector along the x\\ direction. Note that 
at the mean- field level we also obtain explicit expressions for the diffusion matrix and the bias ||23|| : F^ = 1, 
F|| = {1 + e-l^l)/2 and e = 2tanh(£;/2). Of course, these may be modified by renormalization. 

The continuity equation ^ admits both homogeneous and inhomogeneous independent solutions, associated with 
the disordered and the blocked phases. The former is our focus here. To ease comparison with simulation data, we 
choose equal densities for both species: p^(x, t) = p . This solution describes the steady state at the mean-field level. 

At the next order in the fi-expansion, we find a Fokker-Planck equation for the fluctuating part, x^- For our 
purposes, the equivalent Langevin equation is more transparent. At this order, its deterministic part is linear and the 
(conserved) noise is Gaussian. After defining the "reduced" average density p = {1 — j)p and 5 = (3 — 7)/(l — 7), 
and focusing on the fluctuations about the homogeneous phase, the result is: 

aa"(x, t) = £"^(V)x''(x, t) - Vv'^i^, t) , (11) 

where the drift matrix is given by 

(fcPf^-.. f a-p)^T^^ -il-Sp)eT\\d\\ pVrV + p£F||a|| \ , . 

^ ^ " \ pVrV -peF||a|| (l-p)VrV + (l-(5/5)eF||5|| j ^ ' 



and summation over repeated indices is implied in ( pT| ) and in the following. The r7^(x, t) are Gaussian white noise 
terms, with average and second moment: 

(7y,r(x,i))-o, 

(77r(x, t),yf (x', <')) - 2a,f 5(x - x')<5(t - t') , (13) 

where a, P = +, — ; i, j — 1,2, . . . d. Due to the bias, the noise matrices (c,^^) = cr°'^ are diagonal but not proportional 
to the unit matrix: 

'^"'=1". (14) 




4 



Note that cr"-^ is symmetric and due to charge symmetry, we also have cr++ = cr . Similar to T ^, (T°f is diagonal 
and isotropic in the (d — l)-dimensional subspace, characterized by a number . In the absence of the drive, our 
model reduces to an equilibrium system, so that the fluctuation dissipation theorem (FDT) holds. In our case, this 
guarantees cr"^ ct F, or, more specifically, 

<T++ = [p(i-2p) + 7p']r 

p3| . However, when driven, this proportionality does not hold in generic ranges of 7 and p, in that the diffusion and 
noise matrices are renormalized differently by the drive e, similar to the situation in the driven single species case 
p^ . Finally, we point out that there is a correlation between 77+ and ri~ due to the fact that charge exchange is 
allowed. This effect is captured by the matrix cr^ which is expected to be proportional to 7 and negative definite 
for non-zero drive as well. 



A. Steady-state structure factors 



Eqns. (|ll-14) are linear equations which are easily solved in Fourier space. Introducing the Fourier components 
for the fluctuations 

X^{K^)= [ dtd'^xx^{x,t)e-'^'^'+^''\ (16) 



and similar ones for the noise, so that 

(r7nk,c.))=0 

ivT (k, u;)rj^{k', u')) = 2atf [{2nY+H{\^ + \^)5{u + c^')] (17) 
the solution to ( pH ) is simply: 

x"(k,c.) = (L-i)«/'zk77^(k,c.), (18) 

where 

L"'3(k, uo) = £"^(ik) - iuS'^'^. (19) 

Note that, in k space, (£++,£^^) and (£+^,£^+) are complex conjugate pairs. 

Not surprisingly, (x^(k,ct;)) = 0, consistent with the decomposition (^). The two-point correlations of x*(k,a;) are 
just the dynamic structure factors, defined as 

^"^(k, c.) [{2^Y+H{\^ + \^')5{uj + uj')] EE (x"(k, u)x''(k', io')) . (20) 

Using (|l8|) and (p7|), the two independent S"s follow: 

To compare directly with simulations, we need the steady-state structure factors 

5"'3(k) [(2^)'^5(k + k')] EE (x"(k, t)x^{k', t)) , (22) 

which are easily obtained from ( pT| ) by an integration over uj, using the residue theorem and noting that the two zeros 
of det(i) simply correspond to the two stable eigenvalues of C: 



u... = -^^±J<^em^(^^)\ (23) 
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To ensure that the system is within the hnear stabihty region of the disordered phase, we must have Ima'1.2 > 0. 
Since —^Tr{C) = (1 — p)krk is automaticaUy positive definite, we only require det(£) > for aU k ^ 0. CoUecting, 
we find: 

^ I £— P _ ka+ k Rc{C-C-+} 

^^""-iTr(/:) dct(£) ^-iTr(/:) dct{C) ^ > 

so that, with the help of we finally obtain: 

[l-p) k<T++k {kTkf+^j^e^Vlkl 



5++(k) 



(1 - 2/5) krk (krk)2 + 4m2r||fcjj 
~p k<T+-k (krk)2 + ^j^e^Tlkl 
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Re{5'+-(k)} 



(1 - 2/5) krk (krk)2 + 4m2rj|fc|| 



(1 - 2/5) krk (krk)2 + 4m2r||A;| 
(1 - ~p) ka+-k (krk)2 - (T^e^rffcil 



(1 - 2p) krk (krk)2 + 4TO2r||fc2 



(25) 



Tmf 9+-rkn = ~P^^-i^ + m (k^++k)er||fc|| 
^ ^ (l-p)(l-2p) (krk)2 +4m2rj|fc2 

(l-/5)(l-5/5) + p2 (k*T+-k)er||fc|| 



(1 - p)(l - 2/5) (krk)2 + 4m2r||fc2 ' 
To simplify the notation, we have defined a "mass" (in the field theory sense), m , via 



Its role is to mark the linear stability boundary, which, in the limit eL — > 00, is given precisely by m? — . Otherwise, 
for finite eL, the system does not reach the stability Hmit as long as (eL/27r)2 < (1 — (1 — 7)2p)/(l — 2p){{2 — -f)2p— 1) 



is satisfied Thus, it is sufficient to impose rri^ > 0, i.e., p < 2(2— y) ' keep the system in the homogeneous phase. 

Similar to the driven lattice gas |^,^ and the two-species model studied earlier these structure factors are all 
singular at the origin. The singularity takes the form of a discontinuity, either in the function itself or one of its 
derivatives. In particular, both 5"++ and Re{S'^ } are discontinuous, so that the ratios 



limi;||^oS'++(0,fc||) 1-2/5 (1 - (5/5)2 — jra^ — 



and 



limk,^o^++(kx,0) (1- ,5/5)2 -p2 (1 _ <1 _ ^ ^ 

r± i-p Tj 



II p II 

limfe||^oR.e{5'+"(0,A:||)} 1 - 2p l-6p~ TY 



(27) 



limk^^oRe{5+-(ki,0)} (1 - ^,5)2 - p2 l-p _ w 

r_L P 



(28) 



are in general different from unity. In contrast, Im{5^ (k)} vanishes for k ^ in any direction. Here, discontinuities 
occur in higher derivatives. Unlike in the driven Ising lattice gas, these singularities do not simply originate in the 

generic FDT-breaking relation J— 7^ but also in the specifics of this particular driven system, reflected in the first 

factor on the right hand side of (j2^j2^) . It is a monotonically increasing function of p, reaching (X) at p = 2{2^-y) ■ 
a result, the amplitudes of the discontinuities diverge as the system approaches the stability limit of the homogeneous 
phase. 



6 



B. Equal-time spatial correlations 



The equal-time correlation functions G"^{x) = (x" (x'+x, t)x^{^', t)) are just the Fourier transforms of the structure 
factors, 



(29) 



independent of x' by virtue of translational invariance. To simplify the transforms, we introduce some changes in 



notation. First, we rescale the lengths and momenta xj 



x^/rj, ^ x\\/Ti, ^ r^ki , k 

af3 



r becomes the unit matrix. Further, we let 
structure factors in much more compact form: 

5++(k) 
Re{S'+-(k)} 



a/3 



ko-ik 



II 



r,?fcii so that 



(7|| /Fii . After some algebra, we can recast the 



fc2 
kcr^k 



- (ko-^k) 
(ko-^k) 



fc4 + 4m2fc2 
k^ + Am^kf. 



Im{5'+-(k)} = (kcr^k) 



eTS fell 



fc4 + 4m2fc2 ' 



(30) 



where fc =1 k I and 



a = 



a = 



(1 - bpf 



(l-p)((l-<5p)2-p2) 

(l_Jp)2(l_2p) 
(l-p)2((l-5p)2-p2) 

(1 - 5m - 2~p) 



cr++- 



(l-p)((l-<5p)2-p2) 



a = 



a = 



(l-p)((l-^p)2-p2) 

(1 - P)P 

(l-p)((l-<5p)2-p2) 

(l-<5p)(1^2/5) 

(l-p)((l-(5p)2-p2) 

p^l - 2p) 



(1 - 2p) 



(1 - 2p) 



r + + 



(l-p)((l-^p)2-p2) 



L(l-p)2((l-^p)2-p2) 



(1 - 2p) 

(1 - 2/5) 



<r = 



p(2 - (1 + 6)~p) 
{l-p){l-2p) 



r + + 



(l-/5)(l~,5p) + p2 
(l-p)(l-2/5) 



(31) 



For the following, only two major properties of the cr's, rather than their detailed expressions, are important. First, 
they are all positive definite within the stability limit of the homogeneous phase, and second, they are diagonal but 

gencrically not proportional to the imit matrix, since they are related to the noise matrices. Now, referring the details 
of the integrations to Appendix A, the transforms can be carried out exactly. Writing r = |x|, r_L = |x^| and ry = |a;|| |, 
we define 



r(l-i) 

An 2 



^ , ^ cosh(mx||) /m\ 
Fi(x)= ' ^ 



(2^) 



Kd-2 (mr) 



sinh(TOa;||) /m\^ 



d-2 

(^) ' K^imr), 



(32) 



where T{z) is the Gamma function and K,y{z) is the modified Bessel function. The correlations can then be expressed 
in terms of these three functions: 
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G++(x) 

Gr(x) 



1 



er 



Zm 



(33) 
(34) 

(35) 



where are the parts of even or odd in x\\, corresponding to the transforms of the real and imaginary parts 
of S'^ . The full correlation is, of course, G^ (x) — Gj! (x) + G+~(x), reflecting the symmetry (^ of the system in 
the presence of the field. 

Leaving the detailed asymptotic expansions of ( ^3[^5| ) to Appendix B, we only indicate the main features here. 
Clearly, the discontinuity singularities of the structure factors translate into power-law decays of the correlation 
functions. In particular, the first terms in (|3^) and (^J) exhibit the well-known r~''-decay which is associated 
with the breaking of the FDT in the presence of conservation laws In contrast, a discontinuity singularity in 
higher derivatives is not sufficient to produce the r"'' power, as illustrated by the form of G^^, (|35|). The second 
terms in (33) and (|3^ ) and the only one in (|35| ) are more subtle, since the exponentials of the hyperbolic and Bessel 
functions compete. Summarizing our results as mr oo, we find for r± ^ : 



G++(x) oc 



al-ai 



G+-(x)^G+-(x)cx 
G+-(x)(x 



rl -(d-l) 



(36) 



where the ... represent exponentially decaying parts. Here, all three key components of the characteristic non- 
equilibrium power law decays are displayed, namely, the necessity of FDT violation, CT[|/r|| ^ o'±/^±, the dipole 
amplitude associated with strong anisotropy in the presence of a conservation law, and the r"'' itself. We emphasize 
again that the odd part of the cross-correlations, G+~(x), is purely short-ranged. 

Along the field direction, the behavior of the correlations is more complex. Here, a novel power, ('*+^)/^^ emerges, 
which will dominate over the "FDT- violating" r'"^, for all d > 1. Thus, for = we have: 



G++(0,X||), Gr (0,X||)(x -r,, ^ 

(0,2:11) oc -sgn(£a;||)r 



O [ max ■{ r,, , r,. 



d+l 



O r 



(37) 



In the parentheses, we have indicated the nea;i-leading term in the asymptotic expansion of Kd-2 . Surprisingly, 

in d > 3 even this power is still more relevant than the more familiar r~'^. We also emphasize that in (|37|), all 
proportionality constants are positive. Thus, the explicit factors of (—1) carry information about the structure of 
particle clusters in this model. In particular, the sign of G+~(0, a;||) shows that negative charges prefer to be located 
'downfield', rather than 'upheld', from positive ones, as a precursor of the blocking transition. In conclusion, the 
spatial correlations are dominated by the expected r~'^ power law, except along the field, where a novel r^^ ('^+i)/2 
decay takes over. Similar behavior is found if charge exchange is not allowed (7 — 0) ||l9|. Thus, this new power 
law appears to be a generic feature of driven two-species models, associated with the excluded volume constraint and 
the opposite bias. We should note, however, that it can only be generated in the presence of at least one transverse 
dimension, i.e., in d > 1. 



C. The J — 1 case 



All expressions simplify considerably when we set 7 = 1, yet they still capture the essence of this two species model, 
namely, the non-trivial correlations between opposite charges: 
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Re{^+-(k)}Hk<x-k)^^^^ (38) 
r ^ , , _i_ — sen(e) 2™ fell 

M^-(k)}^(k.-k)^l,i±_/, 

where now 4™^ = (1 — 2p)^e^r||. No instabilities can occur here: even for p = 1/2 where the "mass" w? vanishes, 
homogeneous configurations prevail since the model reduces to a driven one-species model. For generic densities, we 
note that the form of the ++ structure factor is the same as in the one-species model, due to the fact that -|-'s cannot 
distinguish between — 's and holes at the microscopic level. The key question is, of course, whether the cr's (especially 
(T++) are proportional to the unit matrix or not. Unfortunately, in the absence of a renormalization group analysis we 
have to rely on simulations to answer these questions. Based on the results of the previous subsection, it is clear that 
only ^^^(k) could possibly produce the r~'^ power law. However, our simulation results indicate that the internal 
symmetry of the system, at this particular value of 7, restores FDT for either species, i.e., the first equation in (|T^). 
This is entirely consistent with the fact that the microscopic steady state distribution of either species is uniform, 
as mentioned in Section II. Thus, correlations will be short ranged, given by a (5-function for identical species and 
exponential decays for opposite charges, except in the field direction, between opposite species, where: 




G-(0,.„).2e(e.„).r^^^(^) ' +0\^\}. (39) 

Here 0{x) is the step function and ajj^" is always negative. Thus, the novel 7-"^''+^)/^ power law, a key feature of this 
two-species model, survives in the cross correlation, even in this simplified case. 



D. Distribution of structure factors 

So far, we have focused entirely on the averages of density-density operators. In this final section, we will con- 
struct the full probability distributions for these fluctuating quantities, i.e., x (k,t)x (-k,t) ^ Re[x (k,t)x (-k,t)] 
Im[x (k,t)x (-k,t)] ^ fQi}Q.^yj]^g ^jjg method of ref. [^. Representing these operators by s^^, s^ and , we seek 
their marginal distributions, for each k-vector separately: 



P++(s++;k) = {d 



X+{Kt)x+*{k,t) 
V 

J lni[x+{Kt)x-*{k,t)] , 

^ [ V 

Here, we have used %*(— k, t) — x^*(k, i), since the densities x^(r,i) are real. Also, we have normalized by V in 
order to obtain a well-defined thermodynamic limit, noting that (27r)'*(S(k = 0) — V. In principle, these distributions 
can be computed explicitly, by inserting the solution x^(k, t) of the Langevin equation ( p]) into ( ^ ) and averaging 
over the noise, associated with (^3|). However, given that these distributions are universal pO(|, depending only on the 
linearity of the Langevin equation and the Gaussian nature of the noise, rather than on the specific forms of diffusion 
and noise matrices, (^2|) and (|lj), a detailed calculation is not necessary. Instead, we can refer to the distributions for 
a simpler case pO| , namely the model without charge exchange, since their forms will be identical to the ones we are 
seeking here. However, some brief comments are in order, to put the results into perspective. For technical reasons, it 
is simpler to compute the characteristic functions (i.e., Fourier transforms) of (M) first. Denoting these by P"^(f2), 
we find that P++(ri) has a single pole in the lower half fi-plane, so that the inverse transform yields an exponential 
distribution for the non- negative variable s++: 

f 1 „-s++/s++(k) -f „++ > n 
P++is++;k)^\ STTJ^^ Its >U ^^^^ 

^ ^ I if S++ < ' ' 
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Here, 5'"'""'"(k) is just the average structure factor, and we will refer to 1/5'"'"+ (k) as the "decay factor" of the expo- 
nential. In contrast, both (Q) and (f2) exhibit two poles, one (ft-) being on the negative, and one on 
the positive imaginary axis, 

= I (Re[5+-(k)] T v/A + (Rc[5+-(k)])2) (42) 

where A = |S'"'""'"(k)p — (k)P > 0. The inverse transforms also result in exponential distributions, characterized 
however by two distinct decay factors |r2+| and 

^, ,„ , +- \ - (43) 
jf el^+l"^ if s+- < ^ ' 

with N = y/A + {Rc[S+-{k)]y. The distribution of Mx-^Wx-'Wl ^ j p+- ^ follows from P+- by just interchanging 
Re[5'"^-(k)] and Im[5'+-(k)] in Eqs. (||,||). 

To summarize, all three distributions are asymmetric exponentials, with representing the most extreme case. 
Due to this structure, their standard deviations are always greater or equal than the averages, so that fluctuations 
will never be "small" , in the usual sense lEol,^ . 



IV. DISCUSSION 



Finally, let us turn to comparisons with simulation results. Typically, we find that power law tails are much more 
difficult to observe than in the single-species case |2^. Apparently, their amplitudes are rather small, so that the data 
are obscured by either critical singularities or finite size effects, depending on the points in the phase diagram which 
we choose to investigate. Thus, we focus on the structure factors. Using a standard least-square routine, we fitted 
our analytical results (eqs. (^0|) and (^8|) before rescaling) to our simulation data. The fit was done simultaneously 
for the three S"s using the smallest 5 x 11 non-zero k vectors. The agreement is quite good, especially considering 
that the theoretical results are based on a linearized Langevin equation, but we note the following: for the 7 = 0.02 
case (Fig. |l|), despite being in the homogeneous phase, the system was relatively close to the continuous transition 
| p6[ , with m A X 10~^ corresponding to a correlation length ^ ~ 25 in units of the lattice constant. In particular, 
"longitudinal" parameters, such as Fn and the cr||'s, seem to suffer considerable renormalizations here. On the other 
hand, the "transverse" parameters, F^ and a±, appear to obey Eqn. (p^). In that sense, the FDT is satisfied within 
the transverse subspace. To illustrate this feature, we combine (Esl), written for the transverse parameters, with the 
explicit form of the structure factors ( |2^ ) for fc|| = 0. This yields the exact "finite size" amplitudes, completely 
independent of k^: 5++(k^,0) = p{l — p) and (k^,0) = —p^, in perfect agreement with the simulations. In 
the full d-dimcnsional space, however, the FDT is of course violated: As a result of the coarse-graining effect in 
the field direction, we generically found cr||/F|| ^CT±/r± for the rescaled "noise" matrices. In particular we had 
CT||/F|| — 0.833 fT^/F^, predicting the typical FDT-violating power law. 

For 7 = 1.00 (Fig. |2|), S'+"*"(k) is completely flat, as we expected, indicating that cr||"+/F|| = a^'^ /Tj_. Moreover, 

the value of this constant is just p{l — p), again consistent with (|l5|). In contrast, (k) clearly exhibits the structure 
of Eqs. (Q). Here the system is far from transitions (^ ~ 3), so that critical fluctuations are completely avoided. 

Consequently, using m — -^(1 — 2p)|e|F| with the mean-field parameters produces a "mass" closely matching the one 
obtained from the fit. 

Now, we turn to a comparison of the analytical results for the structure factor distributions with the simulations, 
summarized in Figs. ^ and ^ for the two smallest wave vectors, respectively. The control parameters were the same 
as those of Fig. [l[ Again, the agreement between our Gaussian theory and the data is quite impresssive. The 
histograms show simple exponential decay ||2^ , while the '-I — ' histograms clearly represent asymmetric exponential 
distributions. To test the theoretical prediction, namely, that the slopes of the histograms are determined entirely 
by the structure factor averages, we simply measured the latter, i.e., S~^^ , Re5'"' and Im5'"' . We then inserted 
the measured averages into the theoretical relations for the decay factors. Clearly, the 'H — h' case is particularly simple 
since the decay factor is just the inverse of 5'"*""'" itself. For the two 'H — ' distributions, the decay factors |ri:p| , given by 
(^), are considerably less trivial, but the agreement is nevertheless remarkable. Here, renormalizations can obviously 
also occur, but can be absorbed into the effective parameters of the theory, leaving the form of the structure factor 
distributions invariant. Moreover, they are automatically captured by the measured structure factors, so that they do 
not spoil the agreement between data and theory here. However, we must avoid critical fluctuations since these fall 
out of the scope of a linear theory. 
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In summary, using both simulations and analytic techniques, we have examined the structure factors in a simple 
model of biased diffusion of two species. We calculated the corresponding spatial correlations, finding not only the 
expected power law decay, r"'', typical for non-equilibrium steady states of conserved systems in the presence of 
strong anisotropy, but also a novel power, r^'"'''^^^^^ , for correlations along the bias, characteristic for two-species 
models. We also investigated the full distribution functions for the structure factors, being universal asymmetric 
distributions. The general agreement between simulations and a Gaussian field theory is surprisingly good, while 
we await a renormalization group analysis of the continuum theory of the model in order to make more detailed 
comparisons closer to the continuous transition. 
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APPENDIX A: MOMENTUM-SPACE INTEGRALS FOR THE CORRELATION FUNCTIONS 



From Eqn. (^0|) we see that we need three basic type of integrals. Although the first one is well known, we list it 
for completeness: 

Then, if cr is diagonal and isotropic in the d—1 dimensional transverse subspace but not a multiple of the unit matrix, 
it is easy to compute 

27r2 7-"+ 

Next, we will outline a formal way to obtain the other two required momentum integrals. For a more rigorous 
treatment see pj. We define Fi and F2 as follows: 



d'^k 



{2t:Y fc4 -I- Am'^k 



It is then helpful to realize that the integrands, without the exponential factor, are simply the convolutions of two 
functions, i.e.: 

d'^h' 

''\F(k')C(k-k') 



-2rm/s|| f d'^k' 



where 



k^ + Am^kf. J {2tt) 



^F(k')5(k-k'), (A4) 



i2 



C(k) = ^^S(k±_) [6{k\\ + im) + 5{k\\ - im)] (A5) 

S'(k) = i^,5(k_L) [5{k\\ + im) - 5{k\\ - im)] . 

The (^-functions with complex arguments should only be understood in an operational sense. The Fourier inverse 
transforms of these functions are easily found: 
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C(x) = cosh(TOa;||) (A6) 
S'(x) = sinh(TOa;j|) . 

Thus, using the convolution theorem, we trivially get 

i^i(x) = cosh(rna;|| )F(r) 

i^2(x) = sinh(ma;||)F(r) . (A7) 

Note that F{r) is the solution of 

{-V^+m^)F{r) ^ J(x) . (A8) 
Then using some algebra and (A8), we can translate VcrV into differentiation with respect to 



VcrVi^i(x) = -<T^5(x) + cr _l2to9|| {smh{mx\\)F{r)} + {a\\ -a±)df^ {cosh(mx||)F(r)} 

= -(Ti(5(x) +CT_L2m9||F2(x) + (cr|| - ct^) ajj Fi(x) 
V<TVi^2(x) ^ aj_2mdi\ {cosh(mx||)F(r)} + (cr|| - a±)df^ {sinh(ma;|| )i^(r)} (A9) 

= a_L2ma||i^i(x) + (o-|| - a^) af[ i^2(x) . 

These forms are particularly useful when we calculate the corresponding long-distance behavior. 
APPENDIX B: LONG DISTANCE ASYMPTOTIC BEHAVIOR OF THE CORRELATION FUNCTIONS 



To obtain the long-distance behavior for V(TVi?(x), we just have to omit the first term in (A2), which is a 
(5-function: 



r(|) r^^ - {d-l)rf, 



This is the typical "FDT- violating" power law, provided that cr is not a simple multiple of the unit matrix. Otherwise, 
the amplitude of this term would be zero. 

Using the "large 2:" asymptotic expansion of the modified Besssel function ||2^ 

we can obtain the long-distance behavior for VcrV Fi{x) and VcrV i*2(x) as m = const. > and r — > cx). Due to the 
strong anisotropics in these functions, we consider three different scenarios: 
(i.) r|| = , r± CO 



Combining (A9) and the asymptotic form of F{r) we find: 



d-l 

m ' 



(27r)t I V^-L 

1 dF{r) 



x«=o '^ii-o r dr 



(B3) 



2; II =0 



d+1 



(27r)t I V^-L. 

while d\\ ^i(x) and dn ^2(x) are simply zero at a;|| =0, since they are odd functions of a;|| . Thus, finally we have 
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d-1 



d+1 



V^VF2(x)|,^_o =0 (B4) 
(ii.) r|| ^ cx) , 7^ : 

In addition to using the asymptotic form of F{r), we can now also write cosh(TOa;||) ~ ^e™'"'i and 
sinh(TOa;||) ~ sgn(a;|| )^e™''ii . In the foUowing, we wiU keep the second leading power in 1/r in order to simplify 
the discussion of case (iii). We find 



r 



d-1 



((.-3)-(..I)^)(^)-.0(-l,)} 



(B5) 



O 



5||Fi(x)~sgn(x||)a||F2(x) 
a^i^2(x)^sgn(.T||)a-[Fi(x) 

Thus, for 7^ we have in leading order: 

/ TT TnTu \ I T*^ \ d.— \ 

+ o(4r)} (B6) 

V<tVF2(x) ~ sgn(x||) VcrVfi(x) . 
(iii.) r|| ^ 00 , = : 

Note that ( p35D was obtained exploiting only r|| — > 00. Setting r j_ = has two important consequences: since now 

r = r||, the exponential decays cancel and, further, the amplitude of the (l/r)~2- term will vanish. Invoking the 
next-to-leading terms in ( p5| ) yields: 

^ ^'"^^-^ (27r)i I 2m Vm/ l^r,]^ 

V^VF2(x)|,^^o ^ sgn(a;||) V^V Fi(x)|,^^o . (B7) 
which are the desired results. 
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(k)} , (c) Im{S'+"(k)} for an i = 100 system at 7 = 0.02 , 
E = 0.279 and p = 0.175. Structure factors are plotted against the integer m\\ = while m± = is taken as a parameter. 



FIG. 1. Steady state structure factors (a) S'++(k) , (b) Re{S'^ 

' and p = 0.175. Structure factors are plot 
Lines are representing the fitted theoretical curves. 
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FIG. 2. Steady state structure factors (a) S'++(k) , (b) Ro{S'+"(k)} , (c) Irn{S+"(k)} for an L = 100 system at 7 = 1.00 , 

E = 00 and p = 0.25. Structure factors are plotted against the integer m\\ = while mx = ^^^j^ is taken as a parameter. 
Lines are representing the fitted theoretical curves. 
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FIG. 3. Histograms representing the distributions of the k = ^{0, 1) structure factors for (a) ''y , (b) — ^—^ — — and (c) 

Im[n^r,_j ^ ^ ^ ^^^^ ^ ^ ^^^^ ^ ^ ^^^g _ ^ ^^^^ Theoretical distributions (a) P++(s++;k), (b) P+"(s+";k) and (c) 
{s^ ; k) are plotted with solid lines on the same graphs. 
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FIG. 4. Histograms representing the distributions of the k = ^i^, 0) structure factors for (a) 

Im[n+n' 
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Re[n+n: ] 

0^ — — and (c) 



,p . L = 100, 7 = 0.02, S = 0.279 and p = 0.175. Theoretical distributions (a) P++(s++;k), (b) P+"(s+";k) and (c) 

{s^ ; k) are plotted with solid lines on the same graphs. 
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